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ABSTRACT 

We analyse the semileptonic decay D + — ► K^-k + 1 + vi using an effective 
Lagrangian developed previously to describe the decays D — ► Plui and D — > 
Vlv\. Light vector mesons are included in the model which combines the 
heavy quark effective Lagrangian and chiral perturbation theory approach. 
The nonresonant and resonant contributions are compared. With no new 
parameters the model correctly reproduces the measured ratio r nres /rnres+res- 
We also present useful nonresonant decay distributions. Finally, a similar 
model, but with a modified current which satisfies the soft pion theorems at 
the expense of introducing another parameter, is analyzed and the results of 
the models are compared. 



I. INTRODUCTION 



The experimental result for the nonresonant semileptonic decay mode 
D+ -> K-n + l + u l is 

r[D+ ^ r ^' 7rV+ ^ (n0nreS 1 ° nant)] = 0.083 ± 0.029 . (1) 

It is important to understand this result before making predictions for the 
other yet unmeasured D — > PiP^lvi decay modes. Furthermore, this could 
also be helpful in understanding of the D — > Plv\ [@, [3| and D — > Vlv\ data 

SI- 

Two main issues arise in the description of D meson semileptonic de- 
cays. The first potential problem is that the D mesons might not be heavy 
enough for the heavy quark effective theory (HQET) to be very accurate. 
The second possible problem is the application of chiral perturbation theory 
(CHPT) in D —>■ Plv\ and D — > PiP^lvi decays, where the light pseudoscalar 
mesons can have quite large energies || |7|, ||, || . 

In order to investigate the semileptonic D13 decays of D mesons we have 
developed a model [[H]] which accommodates the available experimental data 
using HQET and the CHPT description of the heavy and light meson sec- 
tors, respectively. The experimental data for the semileptonic -D/3 decays 
are unfortunately not good enough at this time to empirically determine 
the q 2 dependence of the form factors. What is known experimentally are 
some branching ratios, based on measuring the relevant form factors at some 
kinematical point and assuming a pole-type behaviour for all the form fac- 
tors. The same assumption is also used in many theoretical calculations, for 
example in || and |TT| . 



In our model [ 10 1 the vertices of the processes considered are assumed 



not to change appreciably from their value at the zero-recoil point, where 



they are predicted in the heavy quark limit. However, in our model JED 
the complete propagators for the heavy mesons are used instead of the usual 
HQET propagators. Assuming that these modified Feynman rules can be 
approximately applied to the entire available q 2 region, one can naturally 
understand why some form factors do have a pole-type behaviour and why 
others are mainly flat, which is also in agreement with the predictions of the 
QCD sum rules analysis ||12||. Moreover, in the region where the heavy meson 
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is nearly on-shell, where HQET is most reliable, the HQET prescription and 



our model |10| almost perfectly overlap, providing a simple and consistent 
picture. By calculating the decay widths of all the measured charmed meson 
Di 3 semileptonic decays we found [|Kj] that our model, which is a simple 
modification of the orthodox HQET, worked well, providing confidence in 
extending it to the decays. 

There are some previous theoretical calculations attempting to describe 
the D — > KttIu semileptonic decays using the heavy quark effective the- 
ory (HQET) || [L3|, [yj. Ref. |L3| includes only the light pseudoscalars, while 
for the understanding of the experimental data one must also include the 
light vector mesons. The authors of ref. |L4j considered resonant and non- 
resonant contributions in the overlapping region and indicated that outside 
the resonant region both contributions could be the same order of magni- 
tude. However, no predictions were made. In the present investigation we 
use our simple and instructive model [P^f| , which was quite successful in de- 
scribing the -D/3 decays, to calculate the nonresonant contribution to the 
D + — > K~7T + l + b>i semileptonic decay. The experimental ratio (|1|) is then 
found to be reproduced without introducing any new parameters, which is 
remarkable, considering the simplicity of the model. 

Since the weak current in the model does not satisfy the soft pion limit 
exactly, we have also investigated a modified current which does have the 
exactly correct soft pion limit for comparison, even though the light mesons 
are not particularly soft in the decay. Of course, the strong Lagrangian 
automatically satisfies the soft pion constraint. 

In Sec. II we briefly summarize the strong Lagrangian describing the 
heavy and light pseudoscalar and vector mesons, based on HQET and chiral 
symmetry. In Sec. Ill the two weak currents we consider are given: the usual 
one, whose form is based on HQET, and another, modified to exactly satisfy 
the soft pion theorems, albeit at the expense of additional parameters. In 
Section IV the expressions describing the Dn semileptonic decays are given 
and used to calculate the resonant and nonresonant decay widths, as well as 
some distributions which might be useful in future analyses of experimental 
data. The results are presented in Section V. Finally, a brief summary and 
a few comments are given in Sec. VI. 



II. THE HQET AND CHPT STRONG LAGRANGIAN 
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Our strong interaction Lagrangian [ lOfl , which incorporates both the heavy 
meson SU{2) spin symmetry Hl5| , fl6H and the SU{3)l x SU(3)r chiral sym- 
metry, spontaneously broken to the diagonal SU(3)v [0, describes the heavy 
and light pseudoscalar and vector mesons. A similar Lagrangian, but with- 
out the light vector octet, was first introduced by Wise ||, Burdman and 
Donoghue ||, and Yan et al. [|TI|. It was then generalized to include the 
light vector mesons in , Jl9| , [|(J . 

The light meson sector of the strong Langangian is 



f2 



+ ^[FMF^m, ( 2 ) 



2g v 
where 



= -{rfd^u ~ud^u ] ) , V^ = -{u ] d^u + ud^u ] ) , (3) 

F^u(p) = d^pv - dvp^ + [p^ p v \ , (4) 
u = exp(m/f) , p^ = i(g v /V2)p^ , (5) 

and II and p are the usual 3x3 Hermitian pseudoscalar and vector matrices. 
/ = 130MeV is the pseudoscalar decay constant and gy = 5.9 is determined 
by the values of the light vector meson masses. 

Both the heavy pseudoscalar and the heavy vector mesons are described 
by the 4x4 matrix 



Ha = lil+MK^ ~ Pais) , (6) 

where a = 1,2,3 is the SU(3)v index of the light flavours. P* and P a an- 
nihilate spin 1 and spin heavy mesons cq a having velocity v, respectively, 
and have mass dimension 3/2 so that the Lagrangian is explicitly mass in- 
dependent in the heavy quark limit m c — > oo. Defining H a = 7°if]7° = 
(-^au7 /i + -Pa75)(l+^)/2; we can write the leading order strong Lagrangian as 



3 



£even = C hght + iTl {H a V^ + V^)H a ) 

+ igTr [HaM^)baH a ] + 0Tr [H b v^ - p") ta # a ] 

+ ^Tr (H b H a H a H b ) , (7) 

where V M in the heavy meson kinetic term makes the derivative covariant 
and also ensures that the kinetic term is chiral invariant, since the heavy 
meson field transforms non-linearly under chiral symmetry SU(3) x SU(3). 
The Lagrangian (0) is the most general even-parity Lagrangian to leading 
order in the heavy quark mass and the chiral symmetry limit. The first of 
the four terms in (0) is the kinetic term for the heavy field H a and is thus 
properly normalized. The second term represents the strong interactions of 
the pseudoscalar meson field with the heavy meson field. The third term 
gives the interactions of the light vector mesons with the heavy field. These 
terms involve two unknown parameters, g and (3, which are not determined 
by symmetry arguments, and must be determined empirically. As we will 
see below, only the parameter g will be relevant in the present investigation. 
Finally, the last term comes from the requirement that the Lagrangian 
(0) reduces to Wise's Lagrangian g in the limit gy — > oo. The vector field 
p M then has no derivatives and can be explicitly integrated out. With this 
requirement the coefficient of the last term is fixed. However, this convention 
is irrelevant for our calculations since the vertex with four heavy fields does 
not appear in any diagrams we need to consider. 

We will also need the odd-parity Lagrangian for the heavy meson sector. 
The lowest order contribution to this Lagrangian is given by 

C odd = iXTi[H a a^F^{p) ab H b }. (8) 

The parameter A is a priori free, but we do know that it is of the order 1/ A x 
with A x being the chiral perturbation theory scale. 

III. THE WEAK LAGRANGIAN 
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The weak Lagrangian for the Cabibbo allowed D meson semileptonic 
decays is given at the quark level by 

cf c=AS=1 = -^rMi - [*y*(i - 7 5 )c] , (9) 

where Gf = 1.17 x 10~ 5 GeV -2 is the Fermi constant and V cs = 0.974 is the 
relevant Kobayashi-Maskawa matrix element. 

Of course, as usual, we have to interpret the quark current in terms of 
meson fields. We will present two different models, (A) and (B), for the weak 
part of the Lagrangian: 

Model (A): In the first model, which is based on the HQET approach, we 
assume that the weak current transforms as (3l, 1r) under chiral SU(3)l X 
SU(3)r and is linear in the heavy meson field. Using HQET one can then 
write the most general weak current contributing to D meson semileptonic 
decays to leading order in 1/M and to the next-to-leading order in the chiral 
expansion as 



where 



and 



- «iJ(p - V)\u ] - a 2 J\v a (p - V) a u ] 

+ [JxVa-JaVx-ie^apJ^latip-V^-asA")^, (10) 



Jx = Tr D hx(l-l 5 )H] (11) 



J = Ti D [(l- l5 )H] . (12) 

The first term in fllCf) , i.e. the one proportional to a (= /oy^mo), is 
O(E ), while the rest is 0(E) 0. In the process D -> Plu 0] one takes 



into account only the first term in (10), which is formally the leading term 



In D — » VZz/ |T^] the terms proportional to ai and «2 must be included 
as well, since the diagrams where they appear are of the same order as the 
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diagrams with the terms proportional to a. Actually, the latter diagram has 
also a D*DV vertex (g), which is 0{E). 



Model (B): In the decays D(p D ) Pi{p x )lvi or D(p D ) Pi(pi)P 2 (p2)^ 
one would expect that the part of the amplitude proportional to po is the 
most important since p 2 D = m 2 D while pf = m 2 -C m 2 D . The procedure 
described in (A) takes into account the leading a term in (|T0|) and neglects 
the higher terms which contain derivatives of the light fields and are thus 
formally next-to-leading order terms. However, what is measured is not 
the amplitude but the matrix elements squared, and because the leptons 
are almost massless the part of the amplitude proportional to (pd — HiPiY 
can not contribute. Writing p^ D = (p D — Y^iViY + {^liPiY one sees that 
the formally large part of the amplitude proportional to p^ D contributes to 
the decay width only through the term (JZiViY- So, unless the coefficients 
0^2,3 in (pTO| ) are found to be numerically negligible compared to a they 
can contribute comparably, even if formally they are next-to-leading order 
terms. However, the term proportional to is of higher order and the 
terms in the last line of fllOD do not contribute at all to D + — > K~tt + 1 + i>i 
decay. Consequently, the formal procedure described in (A) does not take 
into account these possibly important contributions. However, if we assume 
that the lower dimensional coefficients of higher dimensional operators in 
( [IDD are naturally smaller (i.e. suppressed by powers of some large scale) 
then we can continue to use the usual approach (A), as was done in the past 

1, eg. 



Another potential inadequacy of approach (A) is that the soft pion the- 
orems are not satisfied by the weak current. This has its origin in the ab- 
sence of all terms of order 1/mo in ([10|). In our case the soft pion the- 
orem requires that the D decay amplitude vanishes in the limit p n — > 
21]. The only term in the current (|TOD that does not satisfy this constraint 
is the first one, i.e. —ifDy/rriDV\Pu\ since it does not have a derivative 
acting on the pion field. In the HQET this term comes from the term 
(/d/\/^_d)-D( 9 — V)\u' ( , where D = e ~ tm Dv.xp. Qne ^ a k es derivative 
of only the exponent and neglects the rest, which is of higher order in l/mr>. 
However, since — (<9 + V)\u^ = A\u^ and the relevant component of A\ in 
our case is proportional to the derivative of the pion field, it is clear how to 
modify the weak current in (A) to satisfy the soft pion theorems: simply mod- 
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ify ([10]) by replacing (P, P*) with (D, D*) and -im D Dv x with D{d - V) A - 
Keeping only the term relevant for our present purposes explicitly one finds 
this modification of ( JTT3| ) to be 



fo 



:D(d - V) x u ] - 2a 3 DAxU 1 - 2a x D{p - V) x u^ - 



2 D(d - V) x (d - V) a (p - V) a v) - f* D ^KEDlv) 



2ao 



m 



(13) 



D 



which we shall refer to as model (B). Note that this modification of the 
original model (A) based on HQET has come at the expense of the appearence 
of a new parameter, viz. 0:3. As a result, the semileptonic D — > Plv data 
do not determine the parameter g in (|7]) as in model (A) [TO, but only a 
combination of g and a 3 . 



IV. THE FORM FACTORS AND DECAY WIDTHS 



Following [|T3J we write down the general form for the matrix element of 
the weak current: 



< n(p n )K(p K ) |s7 M (l - rf)c\D(p D ) >= ir(p D -p K - p w ) 



+iw + (p K +Pn)fi + iw-(p K -Pn)fi- 2he mPl p a D p K pl ■ 



(14) 



The form factor r does not contribute to the decay width if the lepton 
mass is neglected and we will not consider it further. The following combi- 
nations of the remaining three form factors will be particularly convenient 
below: 



F 3 



Xw+ + 



f 1 ( 711 K ~ m w\ 

m D - s Kw - si v ) cos 9 K + I X 



(3X(s K7T s lu ) 1/2 h . 



(15) 

(16) 
(17) 
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Here 6k is the angle between the kaon three-momentum in the Kir rest frame 
and the direction of the total momentum of the Kir center of mass in the D 
rest frame and 



Sk-k 

X 



{PK + Pn) 2 , 8 lv = {p D ~Pk ~ Pn) 2 , 

-[m% + s 2 Kn + s 2 v - 2m 2 D s K7r - 2m 2 D s lu - 2s KlT s lu 



1/2 



(18) 
(19) 



[s 2 Kn + mj c + m*- 2s Kn m 2 K - 2s K7[ m 2 n - 2m 2 K m 2 ^ 2 . (20) 



The D + meson differential semileptonic decay rate can then be written 



as 



d 3 T 



dsKndsi u d COS 9 k 



G 2 F \V CS \ 2 X(5 > 
(4vr) 5 mf ) 3 X 
[\F 1 \ 2 + s\n 2 6 K {\F 2 \ 2 + \Fz\ 2 )] 



(21) 



where the physical region of phase space is defined by | cos6k\ < 1, < s\ v < 
(mo — ^sk-k) 2 and (ttik + rn n ) 2 < sk-k < rn? D . The form factors Fj in (pl|) 
have both resonant and nonresonant parts, which we separate by defining 



F r + F r 



1,2,3. 



The resonant parts are given by 



(22) 



m D 



C 9v 
Cg V y/mBo:i 



,2 \ a l 
jn D - S K n ~ S lu ) — 



Cg v 2X, 



llTT'Ds* ™Ds*fDs* 

m D s lv - m 2 Dsif 



A 



x 2 ■ 

—a 2 



m 



D 



cos 9 



K 



(23) 
(24) 
(25) 



where 



(7 = 8, 



I its It 



ITIk*^ 'k*(SKtt) 



(3 sktt - mL + im K *T K *(sKn) 



(26) 
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and 



2 9o7T mif* 

It is easy to see that for the resonant parts of the form factors, in the zero 
width approximation (IV* — > 0), one obtains the previous expressions for 
the D + -> K*°lvi decay 



The nonresonant contributioin to the decay rate is 

d 3 r nr c||K s | 2 xp 



v 

dsKndsi u d cos 9k (4ii) 5 m 3 D 3 

x [|F™f + sin 2 MI^Tf + |F 3 nr | 2 ) + 

2Re(F 1 nr F[* + sin 2 9 K (F™F™ + F 3 nr F 3 r *))] . (28) 

The nonresonant form factors will be calculated from the leading order 
Feynman diagrams and given in the next section. Note that this nonresonant 
contribution contains not only the nonresonant amplitude itself, but also the 
interference terms with the resonant contribution. 



V. RESULTS AND DISCUSSION 



For the calculation of the Feynman diagrams we use the strong La- 
grangian described in Sec. II and the weak Lagrangian from Sec. Ill for 
both the currents in models (A) and (B). As briefly summarized in the in- 
troduction, and discussed in detail in |1(J, we use the vertices as given by 
our Lagrangian, assuming that they do not vary appreciably away from the 
maximum recoil point, where the HQET is applicable. However, we use the 
complete heavy meson propagators, instead of the HQET approximation. 
The nonresonant form factors can then be straightforwardly calculated from 
the Feynman diagrams with the result: 



9 fp*m D .^m D m D * P*(pd - Pn) 

W + ~ - \9 —2 1 i 1 



fnU [{pd -Pn) 2 - m 2 m ) m 2 m 
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m. 



w 



g fD*m D>l:y /m D m D 



ct2 cos Ok H — — 

SKtt 

Ptt{PD -Pvr) 



ml 



-Oil 



(29) 



fief* [ijpo -Ptt) 2 -rah* 



■[1 + 



m 



+ 



/k"/tt 

-^9 2 fDs*m Ds ^y/m D m Ds 



fxfnKPD -Pw) 2 - m 2 m ) \{p D -Pk-Pk) 2 - m Ds *\ ' 



(30) 
(31) 



In [|TUJ the parameters A, «i and a 2 were fitted to correctly reproduce 
the D + — > K*° decay. In the same reference it was found that the decay 
mode D° — > K~ fixes the parameters g. Due to the nonlinearity of the 
equations involved there are 8 possible sets of these 4 parameters. The same 
values for the decay constants and masses as in plj will be used here, namely 
fo = fo* = (0.24 ±0.05) GeV, f Ds = f Da * = (0.27 ±0.05) GeV and 
m D = 1.87 GeV, m D * = 2.01 GeV , m Dsif = 2.11 GeV H- In deriving ©- 
(|3TD, and in the following, the approximate relation gy = mx*/\ / fKfn = 5.9 
was used for simplicity and will not affect our conclusions. 

The parameters d + and gL depend on the choice of the weak current. 
In model (A), d + = fo and d_ = 0, while in model (B) d + = gL = fo — 
2y/moa 3 . In the soft pion limit, m n — > and p n — > 0, the constraint 



wT + w™ 



2T[ . Combining 



= 
and 



(32) 



indeed, one can easily 
B), where 



should be satisfied 

verify that (|32|) is satisfied in the soft pion limit only in model 
gL|_ = gL, but not in model (A), which was the main motivation for including 
model (B), also, in this analysis. 

We note (^9|)-(pl|) agree with the rsults of Wise, et al. [0, when the 
differences off-shell between our phenomenological approach of model (A) in 
which d + = fo and gL = and the exact HQET [|13j are taken into account, 
as discussed above. Next we present numerical results for both models (A) 
and (B). 



In model (A) there are no unknown parameters. As we discussed above and 
showed in |I0[ , there are 8 possible sets of parameters g, A, cti and a 2 com- 
patible with all the D i3 decay data. For the decay D + — > K~-n + l + i>i we 
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have calculated the ratio R = r(nres)/[r(res) + r(nres)] for all eight possible 
sets and the results are given in Table |l|. We see from Table [l] that all the 
combinations of the allowed values of the input parameters predict a ratio R 
which is consistent with the experimental value (p]) -R exp = (8.3 ±2.9)%. The 
errors quoted are due to the errors only in the model parameters. Unfortu- 
naletly, due to the large experimental error in R the parameters of the model 
are not restricted further. Indeed, the fact that all the 8 sets of parameters 



that were determined from the Di% data |10| give an acceptable prediction 
for R in this Dn decay is remarkable, even given the large uncertainties. 
Perhaps, this merely indicates that R is not very model dependent, provided 
the model fits all the D13 data. 

Since phase space favors smaller Kir energies the tree approximation, 
even though suspect for large energies, is adequate for our purposes. We 
explicitly checked that the nonresonant contribution does not become large 
at high values of sk-k by calculating the distribution d{\ogY nr ) / dsx-n- To 
illustrate this point we show this distribution in Fig. 1, taking as input the 
mean values of the parameters displayed in the 5 th row in Table [I], which 
predict a value of R near to the central value of R exp . In Fig. 1 three curves 
are shown: the dashed line is the contribution of the nonresonant amplitude 
squared, which omitts the last line in (|28|) , the dotted line is the contribution 
of the resonance- nonresonance interferece terms, only the last line in (|28|), 
while the full line is the sum of both these contributions (^). It is clearly seen 
that the nonresonant contribution to the decay width decreases considerably 
at large sk-k- After integrating the curves in Fig. 1 over the invariant mass 
squared of the Kn system, we found that almost 80% of the entire result (f28|) 
comes from the nonresonant amplitude, while only approximately 20% comes 
from the resonant-nonresonant interference term. Fig. 1 is a prediction of 
our model (A), which can be tested in future experiments. 

There is another distribution that is certainly interesting to compare with 
experiment to test this model. It is the distribution of the charged lepton 
energy E\ in the D + rest frame. From the kinematics 

p 1 ( m 2 D + s lv - s K7T \ 

El= 2^{ 2 +X<x»0i) • ( 33 ) 

where 81 is the angle between the charged lepton momentum in the lepton 
pair center of mass system and the direction of the lepton pair center of mass 
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momentum in the D + rest frame. X is defined above in fll^p. 
In general, the distribution in E\ is defined by 



with |TJ 



dSKirdSiyd COS 



2m y) 



dsK-wdsiud cos 6>a' X 



(34) 



d 3 r 



dsK-ndsiyd COS 0# 



(47r) 5 m^ 2 



x 



x { i [l^i| 2 + 2 sin2 ^(l^| 2 + l^3| 2 )]- 

^[\E 1 \ 2 -^sm 2 9 K (\E 2 \ 2 + \E 3 \ 2 )]cos29 l 
J Re(F 2 *F 3 )sin 2 ^cos^} . 



(35) 



In (34) the integration region is defined as in (^), but with the additional 
constraint | cos Q\ \ < 1 where cos Q\ is given by (|33| ) . The range for E\ is 
< Ei < [mjj - (m K + m n ) 2 ]/(2m D ). 

In Fig. 2 this distribution (dT nr / dE{) /Y nr is given by the solid line. In 
addition, the contribution of the nonresonant amplitude alone is given by the 
dashed line, while the contribution of the interference between the resonant 
and nonresonant amplitude is given by the dotted line. 



Model (B) satisfies the soft pion theorems exactly and is defined by the 
current (|T3|). It has the nonresonant amplitudes (p9[-j3T|), but with the pa- 
rameters d + = gL = fu — 2y/rrToa3. Only one combination of d + and g is 
now determined by the semileptonic decay data. The relevant form 



factor 10 becomes 



fx V 2 

and the D i3 decay width is given by |10 



d + m Dsif ^/m D m Ds 
— + 9JDs* 



m 



(36) 



Ds* 



G 2 F \V CS \ 2 
24tt 3 



(m D -m K ) 2 



m D + m K ~ Q 



2\2 



-i 3/2 



4m 2 



m 



K 



D 



(37) 
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After determining the relation between the parameters g and d+ from the 
Diz decay data we then fit the remaining parameter from the non-resonant 
decay data. There are, however, 8 possibilities: four possible sets of 
parameters A, a,\ and a 2 for either of the two possible relations between g 
and d + . From (|36D and (]37f) one can see there is a relation of the form 

-bd + ± ^d\ - 4c« - r) 
9±fDs* = — • (38) 

Labelling the first four sets as given by the first, second, third and fourth 
column in Table [l] with g = g_ by 1 to 4 and the same sets for A, a% and 
ct 2 but with g = g + by 5 to 8, we present their predictions as functions of 
d+ in Fig. 3. All the other parameters are the same as in model (A). The 
experimental band between(_R exp ) min and (-R e x P )max around the central value 
([ip is too broad to determine the model parameter d + = fa — 2 v /moa 3 . 
We can only conclude that model (B) is also compatible with the available 
experimental data ([[]). Further precision is needed to make more conclusive 
statements and discriminate between models (A) and (B). 



VI. CONCLUSIONS 



We have used the effective model developed in ]TD| to calculate the non- 
resonant contribution to the D + — ► K~tt + 1 + vi semileptonic decay. The result 
agrees with the experimental data, giving further support to the model con- 
sidered. More precise experimental data are, however, needed in order to 
reduce the uncertainties in the parameters and further test the assumptions. 

We have also calculated the distributions in the K — ir center-of-mass 
energy and in the charged lepton energy. Both will be useful in comparing 
with future data will provide more sensitive tests of the model. 

Since our original model flO| does not obey exactly the soft pion theorems, 
a slightly modified version, but involving another parameter, was introduced 
and explored. The predictions of this model are also in agreement with 
the present experimental data. Due to the large experimental errors we are 
unable to distinguish between the two models at present. 

We thank Dave Buchholz for very useful discussions. This work was 
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9 


A [GcV- 1 




ai [GeV 1 / 2 ; 




a 2 [GeV 1 / 2 ; 




R 


[%] 




0.08 ±0.09 


-0.34 ±0.07 


-0.14 ±0.01 


-0.83 ±0.04 


6.4 ±0.6 


0.08 ±0.09 


-0.34 ±0.07 


-0.14 ±0.01 


-0.10 ±0.03 


9.0 ± 1.1 


0.08 ±0.09 


-0.74 ±0.14 


-0.064 ±0.007 


-0.60 ±0.03 


4.7 ±0.5 


0.08 ±0.09 


-0.74 ±0.14 


-0.064 ±0.007 


±0.18 ±0.03 


6.4 ±0 


.8 


-0.90 ±0.19 


-0.34 ±0.07 


-0.14 ±0.01 


-0.83 ±0.04 


8.0 ±5.2 


-0.90 ±0.19 


-0.34 ±0.07 


-0.14 ±0.01 


-0.10 ±0.03 


3.0 ±3.5 


-0.90 ±0.19 


-0.74 ±0.14 


-0.064 ±0.007 


-0.60 ±0.03 


11.3 ±6.2 


-0.90 ±0.19 


-0.74 ±0.14 


-0.064 ±0.007 


±0.18 ±0.03 


5.1 ±4.7 



Table 1: The predictions in model (A) for the ratio R = T nr /(T r + T nr ) in 
the decay D + — > K~tt + 1 + i>i. Results are given for all the values of the input 
parameters determined from the Di$ decay data [10] . The experimental ratio 
is = (8.3 ± 2.9)%. 
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FIGURE CAPTIONS 



FIG. 1: The contribution to the Sk-k distribution for model (A) of the non- 
resonant and interference terms and their sums in Eq. (p8|) . The dashed 
line denotes the term given by the square of the non-resonant amplitude 
(nr)x(nr), the dotted line is the interference of the resonant and non- resonant 
amplitudes (r)x(nr), while the full line gives their sum (nr)x(nr) + (r)x(nr). 

FIG. 2: The charged lepton energy distributions for model (A). The dashed 
line denotes the term given by the square of the non-resonant amplitude 
(nr)x(nr), the dotted line is the interference of the resonant and non- resonant 
amplitudes (r)x(nr), while the full line gives their sum (nr)x(nr) + (r)x(nr). 

FIG. 3: The eight possible predictions in model (B) for the ratio R = 
T nr / (r r + r nr ) as a function of the parameter d + = oL = fo — < 2\J m D& , i 
compared with the experimental central value (-R ex p)o an d the values allowed 
by one standard deviation (-R exp )mm and (-Rexp)max- 
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